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ABSTRACT 

Sterile neutrinos are one of the most promising Warm Dark Matter candidates. By 
considering their radiative- and pion-decay channels, we derive the allowed contribu- 
tion of sterile neutrinos to the X-ray, optical and near- infrared cosmic backgrounds. 
The X-ray background puts a strong constraint on the mass of radiatively decaying 
neutrinos {m,ys ~ 14 keV), whereas the allowed mass range for pion-decay neutrinos 
(for a particle lifetime > 4 x 10^^ s) is 150 < mi,s/MeV < 500. Taking into account 
these constraints, we find that sterile neutrinos do not significantly contribute to the 
optical and near-infrared background. We further consider the impact of sterile neutri- 
nos on reionization. We find that the Thomson optical depth due to sterile neutrinos 
is Te = (0.4 — 3) X 10~^ in the case of radiative decays, and it is « 10~^ for the 
pion-decay channel. We conclude that these particles must have played only a minor 
role in cosmic reionization history. 
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1 INTRODUCTION 

Different dark matter candidates can be distinguished on the 
basis of tlieir velocity dispersion, corresponding to a free- 
streaming length, below which dark matter fluctuations are 
suppressed. In Cold Dark Matter (CDM) models particles 
have negligible free-streaming length with respect to scales 
of cosmological interest. On the contrary, in Warm Dark 
Matter (WDM) scenarios the velocity dispersion of the par- 
ticles is sufficient to smear out the fluctuations up to galactic 
scales, depending on the mass of the particles (Padmanab- 
han 1995). For this reason WDM models can alleviate the so 
called substructure crisis, which represents one of the most 
serious problems of CDM theories (Bode, Ostriker & Turok 
2001; Ostriker & Steinhardt 2003). 

Sterile neutrinos are considered the most promising WDM 
candidates (Colombi, Dodelson & Widrow 1996; Sommer- 
Larsen & Dolgov 2001). They can exist only if neutrinos 
have non-zero mass and mixing angles, as predicted by the 
standard oscillation theory (Dolgov & Hansen 2002; see Dol- 
gov 2002 for a complete review of sterile neutrino proper- 
ties). They are massive; so they can decay into lighter par- 
ticles, following a large number of different possible chan- 
nels (Dolgov 2002). Among all the possible decay processes, 
those which lead to the emission of photons are particu- 
larly interesting, because we can put constraints on the 
existence, the mass, the lifetime, and other properties of 
sterile neutrinos, depending on the detection of such radi- 



ation (De Riijula & Clashow 1980; Stecker 1980; Drees & 
Wright 2000; Abazajian, Fuller & Patel 2001a; Abazajian, 
Fuller & Tucker 2001b; Dolgov 2002 and reference therein). 
In principle, radiation due to neutrino decays should be 
seen both as line emission in the local universe (at energy 
Ej — (m^s — rnia) /'^'^vs, where m^s and m^a are respec- 
tively the mass of the sterile and of the light neutrino), and 
as background radiation, given by the integrated contribu- 
tion of decays occurring at different redshifts. Until now, 
none of the experiments built to detect emission lines in the 
local universe has provided any evidence of the existence 
of sterile neutrinos (Henry & Feldman 1981; Shipman & 
Cowsik 1981; Holdberg & Barger 1985; Bowyer et al. 1999). 
The absence of detection can be used to fix upper limits 
on the mass and lifetime of sterile neutrinos. In particu- 
lar, measurements of X-ray emission from galaxy clusters 
(Abazajian, Fuller & Tucker 2001b) provide a strong upper 
limit of the rest mass of radiatively decaying sterile neu- 
trinos (m^3 ^ 5 keV in the zero lepton production mode). 
This limit is particularly constraining, given that studies 
(based on Lya forest data) of power spectrum fluctuations 
due to free-streaming effects provide a lower limit for the 
mass of sterile neutrinos mi,^ 2 keV (Viel et al. 2005; see 
also Narayanan et al. 2000 ^). Studies of cosmological con- 
tributions of sterile neutrinos to the background radiation in 



The lower limit of 0.75 keV for thermal neutrino mass derived 
by Narayanan et al. (2000) can be converted into a lower limit 
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different bands (Stecker 1980; Kimble, Bowycr & Jakobsen 
1981; Ressel & Turner 1990; Davidsen et al. 1991; Abaza- 
jian, Fuller & Tucker 2001b) can only provide upper limits 
of their mass and lifetime. In addition, upper limits derived 
from background radiation can be biased by dark matter 
clustering (Abazajian, Fuller & Tucker 2001b). 

It has also been proposed that radiatively decaying ster- 
ile neutrinos can be responsible of the reionization (Sciama 
1982; Rephaeli & Szalay 1981; Sciama 1990; Sciama 1993). 
This idea docs not encounter particular success today, bo- 
cause there are no observational evidences supporting it. 
The main failure of this model is that the ionizing sources 
proposed by Sciama, sterile neutrinos of mass m^s — 27.4 ± 
0.2 eV, must bo excluded on the basis of the lower limit de- 
rived by the power spectrum analysis (Narayanan et al. 2000; 
Viel et al. 2005). On the other hand, even if it seems unlikely 
that sterile neutrinos are the only responsible of the reion- 
ization, we cannot exclude that they partially contribute to 
it, together with other sources. Furthermore, the high elec- 
tron scattering optical depth (re = 0.17 ±0.04) resulting by 
WMAP data (Spergel et al 2003; Kogut et al. 2003), com- 
bined with the indication that the universe was not fully 
ionized until rcdsliift 6-7, as one can derive from the HI ab- 
sorption present in the spectra of highest redsliift quasars 
in the Sloan Digital Sky Survey (Becker et al. 2001; Fan 
et al. 2003), suggests a long (and, maybe, complex) reion- 
ization history, difficult to achieve considering only "tradi- 
tional sources" (as Population II stars and quasars). One 
of the proposed scenarios (Cen 2003; Furlanetto & Loeb 
2004) predicts a first high redshift phase of partial ioniza- 
tion, triggered by quite exotic objects, followed by a most re- 
cent phase of complete and definitive re-ionization, ended at 
redshift about 6 and due to normal Population II stars and 
quEisars. Among the various candidates proposed as sources 
of the first partial ionization, like massive metal-free stars 
(Haiman & Holder 2003; Choudhury & Ferrara 2004) and 
micro-quasars (Madau et al. 2004), one could reasonably in- 
clude also radiatively decaying neutrinos. 
As an alternative, Hansen & Haiman (2004, hereafter HH) 
have recently suggested that also sterile neutrinos which do 
not directly decay into photons can be a source of this partial 
ionization. They noticed that very massive sterile neutrinos 
(140 < m^s < 500 MeV) mainly decay into pions and lep- 
tons (electrons or positrons) . These electrons are sufficiently 
energetic to Compton-scatter CMB photons, which can then 
ionize hydrogen atoms directly or through secondary elec- 
trons. 

In this paper we want to check the possible contribu- 
tion of sterile neutrinos to the background radiation, in the 
light of the most recent measurements of the extragalac- 

tic background light (EBL), with particular care for the 
optical- infrared (Madau & Pozzetti 2000; Bernstein et al. 
2002; Matsumoto et al. 2000; Wright 2000) and for the X- 
ray band (Gruber 1992; Hornschmeier et al. 2001; Tozzi et 
al. 2001; Moretti et al. 2003 [M03]; Dijkstra et al. 2004 [D04]; 
Bauer et al. 2004 [B04]). We, further, discuss the role of ra- 
diatively and pion-decaying sterile neutrinos in the cosmic 



reionization, calculating the corresponding Thomson optical 
depth. 



2 LIGHT FROM NEUTRINO DECAYS 

In this Section we calculate the possible contribution of ster- 
ile neutrinos decaying at different redshifts to the observed 
background fiux. 



2.1 Radiative decay channel 

To start, let us consider the case of radiatively decaying 
neutrinos. Radiative decay is one of the possible outcomes 
of the following process. 



Va + l+l, 



(1) 



where i>s is the sterile neutrino, Va is an ordinary active 
neutrino and I is a lepton. In particular, this decay is radia- 
tive when the two leptons in the previous expression are an 
electron-positron pair, their annihilation producing a pho- 
ton. This photon has an energy = {mis — rnta) /2mva- 
Stecker (1980) calculated the fiux emitted by sterile neutri- 
nos radiatively decaying at different redshift. Correcting its 
calculations and updating them to present values of cosmo- 
logical parameters^, we find a photon fiux (units of cm~^ 
s-i sr-i): 



I{Eobs) 



Us C 



e-*(^'/"So(5((l + ^)£o6. - So) 



1/2 



d«,(2) 



where ria and r are the present number density and the 

lifetime of sterile neutrinos respectively; Eo and Eobs arc 
the emitted and the observed energy of the photon, and t{z) 
is the time elapsed from the Big Bang to redshift z, which 
can be approximated at high redshift as 

tiz)c^lH^'Q-l/'{l + z)-'^' (3) 

Since S{{l+z)Eobs-Eo) / Offand only if (1 + z) = Eo/Eobs, 
the eq. (2) becomes (for comparison, see Masso & Toldra 
1999) 
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(4) 



This equation depends on two fundamental, substantially 
unknown, parameters; the density and the lifetime t of 
sterile neutrinos. A reasonable upper limit of ris can be ob- 
tained by imposing that all the dark matter is composed by 
sterile neutrinos (Dolgov & Hansen 2002): 



= 1.2 X 



Km^sJ V 0.23 / V0.72y 



(5) 



where Ha is the present density of active neutrinos and ruvs 
is the mass of sterile neutrinos. 

The lifetime r can be expressed as a function of Us and 
Tfii/s , as follows. The rate F for the process (1) is given by 



^ We adopt the following cosmological parameters: Hubble con- 
stant Ho=72 km s-i Mpc"!, Hqm = floM + ^b=0-27, Oa=0.73, 

of ~ 2 keV for sterile neutrino mass, as shown by Hansen et which are in agreement with the recent WMAP determination 

al. (2002). (Spergel et al. 2003). 
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(Boehm & Vogel 1987; Drees & Wright 2000; Abazajian, 
Puller & Patel 2001a): 



5 



(6) 



768 7r3 " 192 TfS 

where 6 is the mixing angle and Gf the Fermi constant. 
The branching ratio for radiative decay, that is the fraction 
of processes of type (1) leading to the emission of a photon, 
can be expressed as (Drees & Wright 2000): 



B 



where Oem is the fine structure constant. Then, the rate 
Frad = "J""^ of radiative neutrino decays is the following. 

9a em 



■ sin^i9 G% m'i 



(8) 



512 7r4 

where sva.'^O can be written as (Dolgov & Hansen 2002^): 

sin^. = 2.8X10- 4^ f ^) '(9) 

where C2 is a numerical coefficient depending on the neutrino 
flavor (c2 = 0.61 for Pe and C2 = 0.17 for p-r and p^) and 
g is the number of rclativistic degrees of freedom at the 
time when sterile neutrinos where produced (wc adopted 
g = 10.75, Viel et al. 2005); finally, gn = sin^^w and gt = 
1/2 + sin^0vK, where sin^^vK = 0.23 (Dolgov & Fukugita 
1992), Ow being the weak mixing angle. 



2.2 Pion decay channel 

We can derive the same quantities calculated for radiative 
decays in the case of pion decay considered by HH. They 

noticed that, if the sterile neutrino has a mass in the range 
140-500 MeV, the dominant decay channel is (Astier et al. 
2001) 



/ + 7r. 



(10) 



where 7r is a pion and I can be an electron, a positron or 
a neutrino. We will focus on the case in which I is an elec- 
tron. It will have an energy _Be = {;inus — m,r) /2, where m,r 
is the pion total mass. Then, because the rest mass of the 
pion is 139.7 MeV, Q < Ee < 180 MeV. As pointed out 
by HH, an electron in this range of energies has a mean free 
path to inverse Compton-scattering with CMB photons gen- 
erally lower than the mean free path to colfisionally ionizing 
hydrogen. Then the most important interaction these elec- 
trons will experience is the inverse Compton-scattering with 
CMB photons. Can these scattered CMB photons signifi- 
cantly contribute to the background radiation? To answer 
this question we made the following calculations. As a first 
approximation, we assume that the electrons produced by 
neutrino decay follow a power law energy distribution with 
spectral index p (this assumption takes into account the fact 
that each electron interacts more than once with CMB pho- 
tons, gradually degrading its energy). The final spectrum 



^ In equation (12) of Dolgov & Hansen 
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the term ( 10^75 ) should be ( 10^75 ) Hansen private 

communication) . 



of photons, initially a Planckian, duo to inverse Compton- 
scattering by a power law distribution of electrons is given 
by (Abramowitz & Stegun 1965; Rybicki & Lightman 1979) 
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where ro is the classical radius of the electron, hpi the Planck 
constant, fcs the Boltzmann constant, T the initial black 
body temperature of photons, p the spectral index of the 
electron distribution and E the final photon energy. F{p) is 
defined by: 



^(j,+3) p^+4p +11 

(p + 3)2(p + 5)(p + l) 



where F is the Euler's Gamma function and ^ is the Rie- 
mann's Zeta function. ^ is a normalization factor defined 
by 



Ay-^A-f 



(12) 



where rie and 7 arc respectively the density and the Lorentz 
factor of the electrons (0 < 7 < 360). 

We included in eq. (11) the dependences on the redshift and 
we integrated it over the line-of-sight, obtaining the follow- 
ing equation. 



I{Eo 



f 

Jo 



dE 



Al 



dVAtAEAQ Az 

.2„2 



Az 



Ho Jo [(1 -h 2)3noM -h Oa] 



1/2 



Az (13) 



where To is the present temperature of the CMB** and 
Eobs = E {1 + z) is the final energy of the photon at red- 
shift z=0. A(z) can be derived from eq. (12). In particular 
we want to express the proper density of electrons ne{z) as 
a function of the initial density of sterile neutrinos ns{z) 
and of their lifetime r. The production rate of electrons is 
approximately given by dne/d£ = ns{z) exp [—t{z)/T]/T. In- 
tegrating, we obtain 



ne(t) 



ns{z) 



fit) At 



(14) 



where f{t) is a function which takes into account the fact 
that the electrons have a finite lifetime. In the simplest case, 
f{t) is a step function 



fit) 



1 if {t(z)-tciz)) <i<t{z) 
otherwise 



where tc{z) is defined as the electron cooling time due to 
repeated Compton-scatterings with CMB photons, tc{z) ~ 
6xl0''(Se/100 MeV)-^[(H-2)/21]-* yr (HH). Substituting 
f{t) in eq. (14), we obtain: 



ne(t) = 



t(z) 



ns{z) 



-i(z)/-, 



At 



{t(z)-tc(z)) 



■ ns{z) e 



1) 



(15) 



We adopted To = 2.275 K (Spergel et al. 2003). 
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where t{z) is computed according to eq. (3). 
Finally ns{z) = ris {1 + z)^ (where Us is the present sterile 
neutrino density) can be easily expressed as a function of the 
present density of active neutrinos ria and of the lifetime r 
and the mass m,^s of sterile neutrinos (HH) : 

A O 

laT"' (16) 



{rrivs/m-KY 



Taking into account that Ua = ^nb/r; (where Ub is the 
present baryon density and 77 is the present baryon-to- 
photon ratio^), we finally obtain: 



A{z) 



4.8(1 + ^)3 Sm 



[{m^s/m^Y - 1] llr? J^""" 7-''d7 



xe-*W/^ (e*^(^)/^ - 1) , 
which can be substituted into eq. (13). 



(17) 



3 BACKGROUND RADIATION 

Now we have all the equations we need to estimate the con- 
tribution of sterile neutrinos to the background radiation. 
Using eqs. (4) and (13) we can derive the observed flux re- 
spectively due to radiative and not-radiative neutrino de- 
cays. In our calculations, we assume that all the photons 
produced (or Compton-scattered) by neutrinos decaying at 
z > 1000 are thermalizcd by CMB photons, and then arc 
not visible. We assume also, for simplification, that all the 
photons produced (or Compton-scattered) by neutrinos de- 
caying at z < 1000 arc not absorbed or scattered. For this 
reason, our calculation represents an upper limit for the op- 
tical and infrared background; however it should give a good 
estimate of the X-ray flux. 



3.1 Soft X-ray background 

A recent estimate of the soft X-ray background (SXRB), 
in the energy range 0.5-2.0 kcV, is provided by M03. Com- 
bining 10 different measurements reported in the literature, 
they found a flux of 2.47 ± 0.11 x 10"® erg s"^ cm"^ sr"\ 
which must be lowered to about 1.48x10"^ erg cm~^ 
sr~^, if the mean contribution of both point and diffuse un- 
resolved sources is subtracted at the level estimated by M03. 
D04 derive the slightly lower value 1.15±1.64 x 10~^ erg s~^ 
cm~^ sr~^, as they subtract a further 1.0-1.7% contribution 
by the diffuse component due to Thomson-scattered point 
source radiation (Soltan 2003). D04 also suggest a maximum 
possible value ~ 4.04 x 10~^ erg s~^ cm~^ sr~^, by sub- 
tracting the lower limit of unresolved sources and of Thom- 
son scattered flux. D04 further notice that the theoretically 
expected amount of X-ray emission in the soft band from 
thermal emission by gas in clusters (Wu & Xue 2001) should 
represent ~ 9% of the total SXRB, considerably higher than 
the 6% estimated by M03. However, D04 do not include in 
their estimate of the background this probable additional 
contribution. Otherwise, the flux due to unaccounted sources 
should be lowered to ^ 4.0 x 10~^° erg s~^ cm~^ sr"'^. 



^ We used n;, = 2.7 x 10"'' cm" 
to WMAP (Spergel et al. 2003). 
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Figure 1. Integrated flux between 0.5 and 2 keV due to radia- 

tively decaying neutrinos as a function of the neutrino mass {solid 
line). The liorizontal dashed line represents the total measured 
background (M03). The shaded area indicates the flux due to 
unresolved sources (B04). 



More recently, B04 investigated the X-ray number counts 
in the Chandra Deep Fields (CDFs). Adopting for the to- 
tal SXRB from 0.5 to 2 kcV the value suggested by M03 
(i.e. 2.47±0.11 x 10"*^ erg s"^ cm"^ sr"^), they found that 
its resolved fraction is 89.515 7%. resolved sources are 
predominantly AGNs (~ 83%) and star forming galaxies 
(~ 3%). This means that unresolved sources produce a flux 
2.59li;4s X 10"^ erg s^^ cm^^ sr^^ 

These strong limits on the SXRB are important to con- 
straint the mass of sterile neutrinos. We base our analysis 
on the paper by B04; in addition we have checked the val- 
ues derived by M03 and by D04. The latter are consistent 
with those reported by B04 and give only slightly different 
constraints on the sterile neutrino mass. In Fig. 1 we report 
the integrated flux between 0.5 and 2 keV due to radiative 
neutrino decays for different neutrino masses. To calculate 
this flux we derived the neutrino comoving density from eq. 
(5), i.e. assuming that all the dark matter is due to sterile 
neutrinos. Given the mass and the density, the lifetime is 
completely determined by eqs. (8) and (9). The flux due to 
sterile neutrinos is compared with the total flux predicted by 
M03 and with the flux due to unresolved sources (B04). As 
we can see from Fig. 1, only neutrinos with masses ^ 950 
keV do not exceed the maximum value of the flux due to un- 
resolved sources (4.00 xl0~^ erg s~^ cm~^ sr~^, B04); the 
same conclusion can be drawn by using D04 data. 

We repeat the same calculation for heavy sterile neu- 
trinos (140 < rrivs < 500 MeV) decaying into pions and 
electrons. In this case we consider the flux due to CMB 
photons Compton-scattered by the electrons coming from 
the neutrino decay (cq. (13), assuming p=l). Fig. 2 shows 
the isocontours of the integrated flux between 0.5 and 2 keV 
due to Compton-scattered CMB photons, as a function of 
the progenitor neutrino mass and lifetime. As we can see 
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Figure 2. Isocontours of integrated X-ray flux in the band 0.5-2 
keV (units of erg cm~^ s"'^ sr~^) due to pion-decaying neutrinos 
as a function of tlic neutrino mass and lifetime. Tiic light gray 
area indicates flux due to unresolved sources (B04). The dark gray 
area indicates the "forbidden" region in which the flux exceeds 
the maximum background flux due to unaccounted sources (B04; 
M03). 

from the eq. (16), the density of neutrinos depends strongly 
on their mass and on the Hfetime. In particular, less mas- 
sive neutrinos have much higher comoving density. For this 
reason, and because of the largely different mass range, the 
observed flux has an opposite behavior with respect to ra- 
diatively decaying neutrinos. In this case the soft X-ray flux 
decreases increasing the neutrino mass. The flux produced 
by 150-500 MeV neutrinos does not violate the SXRB limit 
(B04), if the lifetime is > lO^'^ s. For shorter lifetimes, neutri- 
nos with mass lower than 190 MeV do violate the constraint 
imposed by unresolved sources, according to both M03/D04 
and B04 estimates. 

3.2 Hard X-ray background 

M03 estimated the total background flux in the hard band 
(2-10 keV, HXRB) to be 6.63 ± 0.36 x 10"* erg cm'^ s"^ 
sr~^. They also show that the resolved fraction of this back- 
ground is 88.815^%. Then only the 11.2% of the hard X-ray 
background can be due to unaccounted sources (like sterile 
neutrinos). 

B04 studied the HXRB flux from 2 to 8 keV, using the X-ray 
number counts in the CDFs. They adopted a total HXRB 
flux 5.88±0.36 X 10"* erg cm"^ 8"^ sr"i and found that the 
resolved fraction is 92.6tg;3, dominated by AGNs (~ 95%). 
This means that the unresolved flux is only 4.35l3'8g ^ 10~® 
erg cm~^ s~^ sr"'^. 

We compared the hard X-ray flux due to sterile neu- 
trinos with the background level estimated by B04, and we 
checked for consistency with M03. In particular, for radia- 
tively decaying sterile neutrinos we found that the integrated 
flux in the range 2-8 keV is lower than the background 
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10 100 
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Figure 3. Integrated X-ray flux in the range 2-8 keV due to ra- 

diatively decaying neutrinos as a function of the neutrino mass 
(solid line). The horizontal dashed line indicates the total back- 
ground in the same band; the shaded area indicates the flux due 
to unresolved sources (B04). 

flux due to unresolved sources, as estimated by B04, only if 
rrii^B ^ 14 keV (Fig. 3). If we consider the unresolved flux in 
the range 2-10 kcV estimated by M03, we find a similar con- 
straint for the neutrino mass: m^s ~ 16 keV. This is a strong 
constraint, if combined with the lower limit for the sterile 
neutrino mass ^ 2 keV derived by Viel et al. 2005. We must 
be aware of two uncertainties in our calculations. First, we 
made the assumption that all the dark matter is composed 
by sterile neutrinos. In addition, we have not taken into ac- 
count the clustering of matter at low redshift, which can lead 
to an overestimate of the flux due to neutrinos (Abazajian 
et al. 2001b). 

On the other hand, the integrated flux (in the range 
2-8 keV) of CMB photons scattered by heavy pion-decay 
neutrinos is lower than the upper limit of the background 
due to unresolved sources (B04) for all the neutrino masses 
from 150 to 500 MeV, if the hfetime is longer than ~ 4 x 
10^'^ s (Fig. 4). We obtain the same result by considering 
(in the range 2-10 keV) the flux due to unaccounted sources 
as estimated by M03. For shorter lifetimes, the minimum 
allowed mass increases. For example, if r = 10^^ s, neutrinos 
with masses lower than ~ 215 MeV exceed the upper limit 
of the background due to unresolved sources as derived by 
B04; whereas, using the estimate of M03 , we flnd ~ 200 
MeV. 

In summary, using the HXRB we put much stronger 
constraints on the mass of radiatively decaying neutrinos 
(mi/s ^ 14 keV) than using the SXRB. Instead, in the case 
of heavy sterile neutrinos decaying into pious, all the masses 
from 150 to 500 MeV are allowed, provided a lifetime longer 
than 4 x lO" s. 

To compare our results with those obtained by Abaza- 
jian et al. (2001b), we have repeated this analysis using Gru- 
ber's model of HXRB (1992). We found that, assuming this 
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Figure 4. Isocontours of the integrated X-ray flux in the band 2-8 
keV (units of erg cm~^ sr~^) due to pion-decaying neutrinos 
as a function of the neutrino mass and lifetime. The light gray 
area indicates the flux due to unresolved sources (B04) . The dark 
gray area indicates the "forbidden" in which the flux exceeds the 
maximum background flux due to unresolved sources (B04). 



model of HXRB, the mass of radiatively decaying neutri- 
nos must be m^a ^ 3 keV, in agreement with Abazajian et 
al. (2001b). However, in the following we will use the con- 
straints derived using the estimates by B04 and MQ3, be- 
cause they come from more recent data and do not require 
any assumption about the shape of the HXRB due to unre- 
solved sources. 



3.3 Optical and near-infrared (NIR) background 

A number of factors, such as interplanetary dust scattered 
sunlight (zodiacal light), terrestrial airglow, dust-scattered 
Galactic starlight (diffuse Galactic light), make very diffi- 
cult to obtain a reliable measurement of the NIR and es- 
pecially of the optical background light. Here we briefly 
summarize the most recent measurements of optical and 
NIR background. Deep optical and NIR galaxy counts 
give an estimate of the EBL fraction coming from normal 
galaxies. Madau & Pozzetti (2000) derived the contribu- 
tion of known galaxies in the UBVIJHK bands from the 
Southern Hubble Deep Field imaging survey. In particular 
for the U, V, B and I bands (corresponding to the wave- 
lengths A = 3600, 4500, 6700 and 8100 A) they found a mean 
flux respectively 2.87lEJ-j|, 4.57tg j?, 6.741 J:^! and 8.041 J;^^ 
in units of 10~® erg s~^ cm~^ sr~^. Following a different 
approach, Bernstein et al. (2002) measured the mean flux 
of the optical EBL at 3000, 5500, and 8000 A, using the 
Wide Field Planetary Camera 2 (WFPC2) and the Faint 
Object Spectrograph, both on board the Hubble Space Tele- 
scope, combined with the du Pont 2.5 m Telescope at the Las 
Campanas Observatory. They found for these three band a 



Figure 5. Measurements of the optical and NIR background: 
Bernstein et al. {filled squares); Madau &; Pozzetti [open circles); 
Matsumoto et al. (crosses); DIRBE data reduced adopting the 
Wright model of zodiacal light (filled circles). 



n.etlli in units of 10 



'® erg s ^ cm ^ sr ^, 



mean flux of the EBL respectively 12.0 



+ 17.7 
6.3 



U.dlllt, and 



considerably 

higher than the contribution of the galaxy counts alone. 

The available NIR background data come from the Dif- 
fuse Infrared Background Experiment (DIRBE) on board 
of the Cosmic Background Explorer (COBE) and from 
the Near InfraRed Spectrometer (NIRS) on board of the 
InfraRed Telescope in Space (IRTS). A summary of the 
DIRBE results can be found in Hauser et al. (1998). Mat- 
sumoto et al. (2000) made a preliminary analysis of the 
NIRS data, estimating the NIRB on the basis of the 5 NIRS 
observation days unperturbed by atmospheric, lunar and nu- 
clear radiation effects. Both the DIRBE and the NIRS data 
show an excess in the NIRB with respect to galaxy counts. 
The amount of this excess depends on the model used to 
subtract the contribution of the zodiacal hght from the mea- 
surements. In fact, there are two different models of zodiacal 
light, that described in Kelsall et al. (1998) and that pre- 
sented by Wright 1997 (see also Wright 1998 and Wright & 
Reese 2000). Fig. 5 presents a summary of all the data we 
have considered. 

The measurements of the NIR and, with less evidence, 
of the optical background, seem to indicate the presence of 
an excess with respect to galaxy counts. The most plausible 
explanation of the NIR excess is that the extragalactic back- 
ground light in this wavelength range is due to the redshifted 
UV and optical light emitted by Pop 111 stars (Bond, Carr 
& Hogan 1986; Santos, Bromm & Kamionkowski 2002; Sal- 
vaterra & Ferrara 2003). In particular Salvaterra & Ferrara 
(2003) developed a model of the NIR background which, 
accounting for the most recent predictions of Pop 111 stel- 
lar spectra (Schaerer 2002) and IMF, nebular emission (i.e. 
the radiation coming from the nebula surrounding the star) , 
and Lya photons scattered by the intergalactic medium, is 
able to fit the NIRS data (Matsumoto et al. 2000) and the 
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Figure 6. Background due to radiatively decaying neutrinos with 
mass 14 keV {dashed line), 6 keV (dotted) and 2 keV (solid). 



DIRBE data with both the methods of zodiacal light sub- 
traction. This model is supported also by the analysis of 
the infrared background fluctuations performed by Maglioc- 
chetti et al. (2003). 

In principle, sterile neutrinos could be an additional 
source of optical and infrared light (Stecker 1980; Biller et 
al. 1998). In this Section we want to check this hypothesis, 
calculating the contribution of sterile neutrinos (in the range 
of masses and/or lifetimes allowed by the constraints derived 
from the HXRB) to the optical and infrared background and 
comparing it with observational data. 

In Fig. 6 we plot the contribution to the background 
light by radiatively decaying neutrinos for three different al- 
lowed masses (respectively 2, 6 and 14 keV). We derived 
their density by assuming that all the dark matter is com- 
posed by sterile neutrinos. The lifetime is uniquely deter- 
mined by the neutrino density and mass. We conclude that 
the flux due to radiatively decaying neutrinos in the opti- 
cal and infrared background is completely negligible, being 
at least 8 order of magnitudes lower than the observational 
data (shown in Fig. 5) even in the most favorable hypothesis. 
We notice that the cut-off at low energy is due to the as- 
sumption that photons emitted at 2; > 1000 are completely 
thermalized. The high energy cutoff simply follows from the 
fact that it has to be 2; > 0. 

Fig. 7 shows that also the optical and infrared back- 
ground due to CMB photons Compton-scattered by elec- 
trons which are the product of the decay of massive sterile 
neutrinos^ is negligible with respect to the observed back- 
ground, plotted in Fig. 5. 

From our analysis we have to conclude that sterile neu- 
trinos cannot significantly contribute to optical and infrared 
background. Discarding the hypothesis of sterile neutrinos 



° We assumed spectral index of the electron distribution p=l, 
which maximizes the optical and infrared flux. 



Figure 7. Background due to CMB photons Compton-scattered 
by electrons produced by sterile neutrinos of masses 215 MeV 
(dashed line), 300 MeV (dotted) and 500 MeV (solid). Neutrino 
lifetime is t = 10^^ s; spectral index of electrons p=l. 



on the basis of our results, the only models which can ex- 
plain the excess observed in the NIR background are those 
involving the contribution of Population III stars (Salvaterra 
& Ferrara 2003; Madau & Silk 2005) and/or of high redshift 
miniquasars (Madau et al. 2004; Cooray & Yoshida 2004). 



4 IMPACT ON REIONIZATION 

As briefly discussed in the Introduction, the most recent 
measurements of the Thomson optical depth (Springel et al. 
2003; Kogut et al. 2003; Fan et al. 2003) suggest a complex 
reionization history, difficult to model without invoking the 
contribution of "exotic" sources (like very massive Popula- 
tion III stars and miniquasars). Sterile neutrinos have been 
also proposed as ionization sources, both in the case of ra- 
diative decays (Sciama 1982), and, more recently, of pion- 
decays (HH). In this Section we reexamine this hypothesis, 
using the calculations presented in Section 2. 



4.1 Radiative decay channel 

Let us consider first the possible role of radiatively decay- 
ing sterile neutrinos. To derive the Thomson optical depth 
provided by this process, we calculate the hydrogen ioniza- 
tion fraction, x, due to photons emitted by neutrino decays. 
Under the hypothesis of ionization equilibrium, this can be 
calculated (Osterbrock 1988) as: 



(1 - x)nH 



Atv- 



dN 



dEdAdt 



asdE = x^n„a{T), 



(18) 



where uh ~ nHO+'nH+ is the total hydrogen number density 
(riHO and n_r/+ being the neutral and ionized hydrogen num- 
ber density respectively) ; x = nH+ / riH is the ionization frac- 
tion; Eth ~ 13.6 eV, dN/di?dAdt is the photon flux (units 
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of cm-2 s-i sr-i erg-^), a(r) = 4.18x 10"" (r/10'*K)-0'^26 
cm"^ is the recombination coefficient; (t_b is the photo- 
ionization cross section of hydrogen atoms. 

After some calculations reported in Appendix A, in the 
case of radiatively decaying sterile neutrinos eq. (18) can be 
expressed as: 



Us (1 + z) 



{l-x{z)) Ta{T)nH{Q) Eo Ho[^om{1 + zf + 0.1.^/^ 

ni5o/{l+z) ^ g-t(Eo/Eobs)/^g-^a(,s{-Bo/i5„6,) 

X / —AEob. (19) 



■Etft/{l+2) 



/ Eg \ fipM U O A 



1/2 



where Tabs is defined by eq. (A3). 

So far, we have neglected the possibility that the ion- 
izing photon produces a secondary electron which is suffi- 
ciently energetic to coUisionally ionize additional hydrogen 
atoms. This is an important effect if if E{z) ^ Eth- To take 
it into account, we modify the cross section by introducing a 
factor [l + 4'{x{z))E{z)/Eth], where (ShuU & van Steenberg 
1985) (j){x) = C {l-x")''; for hydrogen ionization C = 0.3908, 
a = 0.4092 and b = 1.7592. Adding this term in eq. (19) we 
finally obtain: 



x'{z) 



Us (1 + z) 



il-x{z)) Ta{T)nHiO) Eo Ho[noM{l + z)'-^ + ^IaV/^ 



So/(i+z) 



Sth/(l+z) 



X f^E(z) 



1+0(X(^)) 



[sioMii+zy-'+nA] 

E,ts{l + z) 



Eth 



- 1/2 

dEobs (20) 



We solve this equation iteratively, evaluating and 
t{Eo/ Eobs) at Zi-1, where [zi — Zi-\) ^ e, with e suitably 
small (in our calculations e = 0.1 gives a good result). 

Finally, we have to include the relic ionization frac- 
tion remaining after recombination. At high redshift this 
fraction can be higher than that due to sterile neutrinos. 
Therefore, to evaluate 4>{x) and Tabs{Eo/Eobs), we take a 
value x{z) = max[xv(z) , Xrei{z)], where Xv{z) is the ioniza- 
tion fraction due to sterile neutrino decays and Xrei{z) is 
the residual ionization fraction. For x,y{z) we use the value 
found by solving eq. (20). For Xrei{z) we adopt the values 
derived by Tegmark, Silk, Rees Blanchard, Abel & Palla 
(1997) from z =1100 to z—20, and the asymptotic, constant 
value Xrei{z)=3xl0~'^ for z < 20. We found that the dif- 
ferences obtained considering other estimates of Xrei{z) - 
as those derived in Jones & Wyse (1985), Seager, Sasselov 
& Scott (2000), or ShuU (2004) - are negligible. The ion- 
ized fractions we calculated in this way are shown in Fig. 8. 
The ionization fraction due to radiatively decaying neutri- 
nos of masses rrius ^ 14 keV (higher masses are prohibited 
by X-ray measurements) is always x < 0.1. For example, 
the ionization fraction produced by a sterile neutrino of 2 
keV (14 keV) exceeds the relic ionization fraction at redshift 
2 ~ 50 (2 ~ 200) and increases up to s ~ 0.003 {x ~ 0.01) 
at redshift z — 10. 

It is then possible to derive the Thomson optical depth 
Te due to radiatively decaying neutrinos using the standard 
expression 



0.1 r 



T 1 I I I I III 1 1 I I I I III 1 1 I I I I II 




Figure 8. Ionization fraction {top panel) and Thomson optical 
depth [bottom panel) due to radiatively decaying neutrinos of 
masses 2 (solid line), 4 [dotted line), 8 [dot-dashed line) and 14 
keV [dashed line). The solid thin line indicates: in the top panel 
the relic ionization fraction (Tegmark et al. 1997); in the bottom 
panel the Thomson optical depth due to relics. 



Te = / Az — C(JTne[z), (21) 

Jo 

where or is the Thomson cross section, At/dz = 
H^^[^om{1 + zf_+ nA]"^'"^ and ne.[z) = nH+{z) = 
x{z)nH{Ci) [1 + z)'^ . The results are shown in Fig. 8. For 
an integration upper limit z — 1000 (i.e. corresponding to 
the low redshift boundary of the last scattering surface), 
we find the results reported in Table 1. The contribution of 
radiatively decaying neutrinos to ionization fraction is neg- 
ligible, being in the range 3.77 x 10"^ < < 2.86 x 10"^ 
for = 2—14 keV, i.e. much lower than the WMAP best 
value (r, = 0.17 ±0.04). 



4.2 Pion decay channel 

To calculate the ionization fraction due to CMB photons 
Compton-scattered by electrons produced by heavy sterile 
neutrinos we use again eq. (18). In this case: 

1 87rV§ 



AN 
AEAAAt 



4.8 



3n6 



1 



11 



7? J^'""" 7-Pd7 Ho 



(l + 2)e-'W (etc{^)/-_l) 
[noM{l + zy^ + ^aY'^ 

where Taba{z, Zdecay) IS giveu by the eq. (A3). Following 
the same procedure as in previous Section (see also Ap- 
pendix A), we find: 



4tt 



[l-x{z)) a{T)nH{Q)(l + z)^ 



dN 



Eth/(l+z) 



AEdKdt 
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X O-B(z) 



Eobs{l + z) 



Eth 



AEobs, (23) 



where dN/dEdAdf comes from the eq. (22), Q-{T), o"_b{z), 
Eth and (j){x{z)) are the same eis defined for radiatively de- 
caying neutrinos. We adopted hfetime r = lO^^s as sug- 
gested by HH. Because of the constraints derived from the 
HXRB, only masses 215 ^ rrivs ^ 500 MeV are allowed as- 
suming this lifetime. 

The results shown in Fig. 9 imply that the ionization 
fraction due to pion-decaying neutrinos is negligible at high 
redshift, being always lower or close to the relic one for 
z ^ 30. At redshifts lower than 20, the contribution of ster- 
ile neutrinos becomes more important, albeit Xe^O.l; as a 
consequence, the Thomson optical depth (Fig. 9 and Table 
1) is very low (re < 10"^). Even for m^s = 150 MeV with 
a lifetime r = lO'^^ s (which is however forbidden by HXRB 
limits), we only obtain re = 5.8xl0~^ (Table 1). This result 
discards pion-decaying neutrinos as a source of reionization, 
in agreement with the results of Pierpaoli (2004). 



4.3 Comparison with previous results 

The values of the ionization fraction and of the Thomson 
optical depth derived in Section 4.2 partially disagree with 
the scenario drawn by HH, who suggest an early ionization 
due to pion-decaying neutrinos. In this Section we clarify 
the possible sources of such discrepancy. 

Let us consider first the fiux of ionizing photons in the 
two models. This fiux can be derived by integrating eq. (22) 
from Eth/i'i + z) to Eem/il + z). At redshift « = 20 (the 
reference redshift in HH) we find that the ionizing photons 



fiux is ~ 30 cm 



In the case studied by of HH, the 



fiux of ionizing photons can be derived from their equation 
(10), i.e. Ub — exnf"'', where nj, is the number density of 
(re)ionized hydrogen atoms, x = Ee/Eth, n^^"'^ is the num- 
ber density of electrons produced by neutrino decays, and 
e ~ 0.3 is the fraction of the initial energy of the secondary 



Figure 9. Ionization fraction {top panel) and Thomson optical 
depth {bottom panel) due to pion-decaying sterile neutrinos of 
masses 500 {solid line), 300 {dotted line), 215 {dashed line) and 
150 MeV {dot-dashed line). The solid thin line indicates: in the 
top panel the relic ionization fraction (Tegmark et al. 1997); in 
the bottom panel the Thomson optical depth due to relics. 



photoelectrons spent in coUisional ionizations (ShuU & van 
Steenberg 1985). In particular, we can write: 



prod 



(24) 



At z = 20, for Ee = 100 MeV, / ~ 9.7 x 10^ cm" " s 
is a difference of more than 2 orders of magnitude with re- 
spect to our result. However, if we take e = x = 1, we obtain 
7 ~ 8.0 X 10"^, which is more than 3 orders of magnitude 
lower than our value. Then, the assumption that an elec- 
tron of energy E^ is able to produce ~ e Ee/Eth ionizations 
is crucial, and, we think, quite optimistic. For example, a 
CMB photon Compton-scattered by such electron at 2: 20 
will have an energy much higher than Eth, and, because 
as oc {E/Eth)~^, it is unlikely that this photon immedi- 
ately interacts with hydrogen atoms. In addition, this high 
energy photon can be involved in many other processes, like 
pair production, and its energy does not entirely contribute 
to reionization. 

However, even under the assumption that a primary 
electron yields x ionizations, additional physics must be con- 
sidered. HH (j) do not consider the effect of recombinations, 
{a) assume that the density of sterile neutrinos is constant, 
(ill) neglect the gas clumping factor, {iv) do not weight the 
photoionization rate by the cross section (Je{z)- We therefore 
add these processes to HH's calculation. 
The number of recombinations (Nrec) and of ionizations 
{Nion) per unit volume can be written as: 

iV.ec = aCx2ni(l + 2)'^ 

N,on = {l-x)^Us€X ^ ' 

where a — 4.18 x 10"^'' cm"^s~^ is the recombination co- 
efficient, C = {njj) / {nn}'^ i ^ 1) is the clumping factor, x 
the ionization fraction and Ub the present baryon density, x 
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and e are the same as previously defined, fats is the frac- 
tion of the total radiation which is absorbed at each red- 
shift (to take into account the photoionization cross sec- 
tion). Finally, due/dt is the production rate of primary elec- 
trons, i.e. drie/dt = neT~^{l + z)'' e"'^^'''^, whore Ue = 
A.8naT~^ {(rnvs/fnTr)^ — and Ua = 3/11 nt??"^ cm~^ 
(see Sec. 2.2). 

Then, imposing iV^ec = Nion we find: 

^ 1.31 e-"-^^/- fabsex , . 

(1-x) aCmv{i + z)^Tmm.,/m„)^ -1)' ^ ' 

where t[z) is given by eq. (3). For comparison with HH, 
wc adopt T — 10^'"' s and z = 20. At this redshift, we can 
neglect^ {1 — x), because a; ^ 1. As iJe = 1/2 {m^a—rriTr) and 
choosing rrivs =215 MeV (i.e. the minimum mass allowed by 
the HXRB for r = 10^^ s, corresponding to the maximum 
a;), we obtain: 

flilC-"^ (27) 

Even if fabs = 1 and C = 1, an ionization fraction x{z = 
20) ~ 0.07 <C 1 is found, i.e. incomplete reionization. More 

realistically, C will be higher than 1, whereas fahs < 1- For 
example, by calculating the difference between the ioniz- 
ing flux due to Compton scattered CMB photons (see Sec. 
4.2) with and without the term e~'^'^*' (see eq. (A3)), we 
find faba ~ 0.1 in the range 500 > « > 10. Then, we have 

x{m^a = 215 MeV) ~ 0.022 {^f^ {^Y^'' C'^'"^ . This 
value is in good agreement with the results of Pierpaoli 
(2003), but it is a factor 10 higher than that reported in 
Fig. 9. This difference can be explained with the fact that 
equation (26) still contains some simplifying assumptions, 
as the term x (see previous discussion), instead of a more 
realistic spectrum for the ionizing radiation. 



5 SUMMARY 

We have studied the possible contribution of sterile neutri- 
nos decays to the cxtragalactic X-ray and infrared-optical 
background, deriving limits on their masses. Such con- 
straints have been then used to assess the impact of ioniz- 
ing radiation produced by the decays on cosmic reionization. 
Both radiatively decaying neutrinos and heavy neutrinos de- 
caying into pions and electrons (HH) have been considered. 

By calculating the contribution of sterile neutrinos to 
X-ray background, we have put a strong constraint on the 
neutrino mass. In particular, by requiring that the esti- 
mated flux due to sterile neutrinos does not exceed the 
SXRB (Moretti et al. 2003; Dijkstra et al. 2004; Bauer et 
al. 2004), we found that radiatively decaying neutrinos must 
have masses nivs ^ 950 kcV. Much more stringent limits can 
be found from the comparison with the HXRB. In fact, 
by requiring that the flux produced by sterile neutrino de- 
cays does not exceed the HXRB (Bauer et al. 2004), we 
find that the mass of radiatively decaying neutrinos must 
be rriva ^ 14 keV. This constraint on radiatively decaying 
neutrinos is very important if combined with the results by 
Viel et al. (2005), who found a lower limit mva ~ 2 keV 

This results in a diflference of the order of 10~* 



from the study of power spectrum. Pion-decaying neutrino 
masses in the range 150 MeV rriva ^ 500 MeV are allowed 
if the lifetime is longer than 4 x lO^'^ s. For shorter lifetimes, 
the minimum mass increases: for example, if r = 10^^ s, 
215 < m^a ~ 500 MeV. 

We have calculated the sterile neutrino decay contribu- 
tion to the optical and infrared background light. We find 
that both radiatively and pion-decaying neutrinos give a 
flux in the optical and NIR range which is several orders 
of magnitude lower than the observed one. Then alternative 
sources, as very massive Pop III stars (Salvaterra & Ferrara 
2003), must be invoked to explain the detected NIR back- 
ground excess. 

Decaying sterile neutrinos might also be a potential 
source of cosmological reionization. We derived the ion- 
ization fraction and the Thomson optical depth produced 
by radiatively decaying sterile neutrinos and CMB photons 
Compton-scattcrcd by electrons produced by heavy (pion- 
decaying) neutrinos. Radiatively decaying neutrinos produce 
an optical depth Te = (0.4 — 3) x 10~^, whereas the con- 
tribution to reionization of pion-decaying neutrinos is even 
smaller, yielding an optical depth Te « 10~^. In conclusion, 
our calculations suggest that both radiatively- and pion- 
decaying neutrinos are not viable reionization sources. 
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APPENDIX A: IONIZATION FRACTION BY RADIATIVELY DECAYING STERILE NEUTRINOS 



The hydrogen ionization fraction, x, due to photons emitted by neutrino decays can be derived, under the hypothesis of 
ionization equiUbrium (Osterbrock 1988), as: 



(1 — x) uh / 47r 



poo 



dN 



d£;dAdi 



ue dE = x'^ n^a(T), 



(Al) 



where TiH = riHa + nH+ is the total hydrogen number density (nno and nH+ being the neutral and ionized hydrogen number 
density respectively); x = nu+jnu is the ionization fraction; Eih = 13.6 eV, dN/dBdAdt is the photon flux (units of cm~^ 
s^^ sr"^ erg-^), a(T) = 4.18 x 10-^^(r/10*K)-°-^2® cm^ s ^ is the recombination coefficient; ob is the photo- ionization cross 
section of hydrogen atoms. 

We have derived dN/dEdAdt from eq. (2). The specific flux will be: 



dN 



1 c na(z) 
cLEdAdt ~ '^H{z) r 



-t(z')/T - 



d{{l + z')E{z)-Eo) 



{1 + Z') [(l + «')^«MW+nA]l/2 



dz', 



(A2) 



where, in the reference frame of the hydrogen atom, z' and are respectively the redshift at which the sterile neutrino decays 
emitting a photon and the redshift at which this photon ionizes the hydrogen atom. Passing from the reference frame of the 
hydrogen atom to the reference frame of an observer (at redshift 0), the rcdshifts z' and respectively correspond to Zem and 
z, through the relation (1 + z') = (1 + Zem)/{'^ + z). As in Section 2, Eq is the energy of the photon when it was emitted. 
Tahs{z') takes into account the fact that a photon emitted by a sterile neutrino will be absorbed, due to hydrogen photo- 
ionization, after a certain mean free path \h- It is given by: 



rabs{z')= I \H{z)^dz 
Jo 



(A3) 



where dZ is the proper distance element. The mean free path of the photon, \h, can be expressed as: 

" nHo{z)aE(z) UHoiO) + z)^{l - x{z))aE(z)' 

Because 71^(2) = (1 + z)^ and taking into account that 5((1 + z')E{z) — Eq) 7^ if and only if (1 + z') = Eo/E{z), eq. (A2) 
becomes 

1 c ns{l + zf 



dN 



dEdAdt AnHiz) 



{^y^M{z)+QA 



1/2 



(A4) 



Substituting E{z) — Eobs{l + z) {Eobs being the energy of the photon today), H{z) = i/o[fioM(l + z)^ + Oa]'^'^^ and 
^m{z) = 'Itl^^ we obtain 



[f20iM(l + 2)^+"A] ' 

dN 1 



Us (1 + zfe 



dEdAdt At: Ho[Q.om{1 + zY + Q.kYI'^ 



tEq 



( En \ (ipM I o 



-1 1/2 ' 



(A5) 



where, assuming that Q.a << {1 + z')^Qm{z), Oa < {1 + z)^floM and ilM{z) ~ 1, and substituting (1 + ^') = {1 + Zem) / 0- + z) 

and (1 + Zem) = [Eo/Eobs), we can find 
2 . 



-3/2 



t{z') ~ -Ho'Q-^/\Eo/Eot 
Similarly, we can also show that 



: t{Eo/Eobs) 



abs{z') = / Xh^{z) —dz = / Ah^(z) — dz = rabs(Eo/Eobs) 
Jo ^ J z ^ 



Finally we can write 

dN 



47r- 



dEdAdt 



as dE - 



c 

Ifo 



Eq 



Us (1 + zf 



CE(z) e"*'^"^^"*"^'^^ g-'rabs(Eo/Eoba) 



T[noM(l + 2)3+nA]V2 



Eo 



Eg 



^^i)M L ( ) , 



1/2 



dE{z) 



(A6) 



(A7) 



(A8) 



where the upper limit of integration, Eq, takes into account that the photon energy at redshift z cannot be larger than at 
redshift Zem, when it was emitted. We can make a change of integration variable {E{z) = (1 + z)Eobs), obtaining: 



Ait 



dN 



dEdAdt 



as dE 



hJe 



EQ/(i+z) 



71,(1+2)" 



,/(l+.) r[fi0M(l+2)3+nA] 



1/2 



Eo 



( Eq \ ^OM 
Eobs J l^OM(i + ^y 



+ f2A] 



1/2 



dEob 



(A9) 
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For c^E(z) we used the following approximation (Ostcrbrock 1988): 



13 



Eth 



+ (1-/?) 



Eobs{l + z) 



-("+1) 



(AlO) 



where ath = 6.30 x 10"^^ cm^ /3=1.34 and s=2.99. 

Finally, substituting eq. (A9) into eq. (Al), reordering, and substituting nniz) = n/f(0)(l + «)^, we obtain the expression 
for the ionization fraction: 



Us (1 + Z) 



{l-x{z)) ra(r)nw(0) Eq Ho [QoMil + z)^ + ^aY/'^ 



( Eq \ ^OM L Oa 



(All) 



